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ABSTRACT 


Cygnus X-1 is a well-studied persistent black hole X-ray binary. Recently, the three parameters needed to 
estimate the black hole spin of this system, namely the black hole mass M, the orbital inclination i and the 
source distance D, have been updated. In this work we redetermine the spin parameter using the continuum- 
fitting technique for those updated parameter values. Based on the assumption that the spin axis of the black 
hole is aligned with the orbital plane, we fit the thermal disk component to a fully relativistic thin accretion 
disk model. The error in the spin estimate arising from the combined observational uncertainties is obtained via 
Monte Carlo (MC) simulations. We demonstrate that, without considering the counteracting torque effect, the 
new spin parameter is constrained to be a. > 0.9985 (36), which confirms that the spin of the black hole in 
Cygnus X-1 is extreme. 


Keywords: X-rays:binaries — black hole physics — binaries:individual (Cygnus X-1) 


INTRODUCTION 


The spin parameter is a fundamental physical property of a black hole, which can be estimated observationally using one of 
two widely-used methods; the iron-line fitting method (Fabian et al. 1989, Reynolds & Nowak 2003) and the continuum-fitting 
method (Zhang et al. 1997). In this work we use the continuum-fitting method to estimate the dimensionless spin parameter ax 
(a, = a/M = cJ/GM?, where M and J indicate the BH mass and angular momentum, respectively) of the stellar-mass black 
hole Cygnus X-1. Up to now, this method has been extensively employed in roughly a dozen stellar-mass black holes: 4U 1543- 
47 (a, — 0431037, 90% confidence, Morningstar & Miller 2014), A0620-00 (a, = 0.12 + 0.19, lo, Gou et al. 2010), Cygnus 
X-1 (a, > 0.983, 3c, Gou et al. 2014), GRO J1655-40 (a, = 0.70 + 0.10, 1c, Shafee et al. 2006), GRS 1915+105 (a, > 0.95, 
3c, McClintock et al. 2006), LMC X-1 (ax = 0.92*005, 1e, Gou et al. 2009), LMC X-3 (ax = 0.25* 0 15, 1a, Steiner et al. 
2014), M33 X-7 (ay = 0.84 +0.05, 1c, Liu et al. 2010), Nova Muscae 1991 (a, = 0.634916, lø, Chen et al. 2016), XTE J1550- 
564 (ax = 0.34*025, 1c, Steiner et al. 2011), IC10 X-1 (a, Z 0.8, 90% confidence, Steiner et al. 2016), GX 339-4 (a, < 0.9, 
Kolehmainen & Done 2010), H1743-322 (a, = 0.2 + 0.3, 1c, Steiner et al. 2012). 

In applying this method, the key fitting parameter is the inner radius of the thin accretion disk Rin, which can be determined 
by fitting the X-ray continuum spectrum to the Novikov-Thorne thin disk model (Novikov & Thorne 1973). Riy is considered 
as equivalent to the radius of the innermost stable circular orbit Fisco (Shafee et al. 2008, Penna et al. 2010, Noble et al. 2010, 
Kulkarni et al. 2011). In addition, Fisco is a monotonic function of the dimensionless spin parameter a,, decreasing from 6 
GM/ c? to 1 GM / c? as the spin increases from a, = 0 to a, = 1. Therefore, given the inner radius, one can constrain the value 
of a, through this simple monotonic relationship (Bardeen et al. 1972). 

In order to make sure the inner radius of the disk has reached the ISCO, the high/soft state spectra dominated by the ac- 
cretion disk component are usually rigorously selected to avoid confusion introduced by the strongly Comptonized component 
(Shafee et al. 2006, McClintock et al. 2006, Liu et al. 2008). The method requires a thin disk with a moderate luminosity, which 
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means that the ratio of the disk thickness and the local disk radius should satisfy H/R < 0.03, and the bolometric Eddington- 
scaled luminosity | = L(a,, M)/ Lgaa < 0.3 (McClintock et al. 2006, Zhu et al. 2012, McClintock et al. 2014). However, in 
practice, the high/soft spectra are found to be difficult to select for some sources, such as Cygnus X-1, which has never reached 
a canonical high/soft state according to its standard definition (Gou et al. 2011). Therefore, Steiner et al. (2009b) proposed an 
alternative approach to select the spectra, namely that the scattering fraction fse < 25%. The scattering fraction refers to the 
fraction of the thermal seed photons that are scattered into the power-law component, and is quantified via the parameter fsc in the 
Comptonization model SIMPL (Steiner et al. 20092). The maximum value of the fraction is assumed to be unity. This approach 
has been successfully applied to several black hole systems (Gou et al. 2009, Steiner et al. 2011, Chen et al. 2016), and it was 
found that some steep power-law (SPL) state spectra could also be used for the spin measurement, in addition to the standard 
high/soft state spectra. This criterion thereby allowed the spectral sample to be extended beyond the pure thermal state. The 
advantage of this approach is not only that all the selected spectra give the roughly same inner radius (i.e., the spin parameters 
remain quite stable both in the thermal state and the more strongly Comptonized SPL state), but also that the approach can be 
readily implemented (Gou et al. 2009, Steiner et al. 2010, Chen et al. 2016). 

The success of the continuum-fitting method relies on having accurate values of three input dynamical parameters beforehand, 
namely the black hole mass M, the disk inclination ¿ and the source distance D. The spin error introduced by the combined 
observational uncertainties of M, 7 and D, which together dominate the error budget in spin, is typically estimated via Monte 
Carlo simulations (Liu et al. 2008, Gou et al. 2011, McClintock et al. 2014). 

Cygnus X-1 is the first confirmed black hole candidate system, and is one of the most famous persistent black hole X-ray 
binaries, whose dynamical parameters have been well-estimated. A systematic study of the system was carried out by Reid et al. 
(2011), Orosz et al. (2011), and Gou et al. (2014). Reid et al. (2011) estimated the distance using trigonometric parallax, pro- 
viding a precise, model-independent distance measurement for the system. Afterwards, using the distance as one of the con- 
straints, Orosz et al. (2011) fit the optical photometric light curves and radial velocity curves to determine the mass of the black 
hole and the inclination angle of the system. With the system parameters in hand, Gou et al. (2011) firstly found a near-extreme 
value of the spin, a, > 0.95 (30), via the continuum-fitting method, and Gou et al. (2014) later confirmed this extreme result but 
with a more stringent lower limit a, > 0.983 (3c). 

Recently, Miller-Jones et al. (2020) made a new set of six daily orbital-phase resolved VLBA observations, which revealed 
an orbital phase dependence in the measured astrometric positions. Hence, free-free absorption in the stellar wind will affect 
the astrometry by changing the position of the optical depth unity surface seen as the radio core, which if unaccounted for can 
affect the measured parallax signal. They therefore combined the archival data of Reid et al. (2011) with new VLBA data, and 
conducted their astrometric fit only in the dimension perpendicular to the jet axis on the plane of the sky, thereby negating the 
impact of the stellar wind absorption. This gave a refined measurement of the distance, of D = 2.22 Es kpc. The updated 
distance measurement required a revision of the optical modelling, leading to new values for both the black hole mass and orbital 
inclination, of M 221.2 + 2.2 Mọ, i = 27.51? ui , with all uncertainties quoted at 1c. Using these up-to-date values, we are 
able to refit the X-ray spectra used in Gou et al. (2011) and Gou et al. (2014), and refine the constraints on the black hole spin. 
The updated results are presented in this paper. 

This paper is organized as follows. In Section 2, we discuss our observations. In Section 3, we describe our detailed spectral 
analysis. A discussion is presented in Section 4. Finally in Section 5, we summarize our results and give a conclusion. 


OBSERVATIONS 


Because Cygnus X-1 has never reached the classically-defined high/soft state, we used the scattering fraction criterion outlined 
in Section 1 to select the spectra. We found several spectra that could be used for the spin measurement, which were originally 
presented by Gou et al. (2011, 2014) (hereafter GOU11 and GOUI4). In this paper, we used exactly the same spectra for our 
analysis. The log of these observations is presented in Table 1 for convenience. 

The 1996 archival ASCA&RXTE spectrum (hereafter SP1) was discussed in GOUI1. We disregarded the other two spectra 
analyzed in GOUI1 because they have fs: above 25% (30.596 + 1.2% and 30.6% + 0.6%, respectively), which is the upper limit 
for reliable measurements using the continuum-fitting method. Details of the 2011 archival Chandra&RXTE spectra (hereafter, 
SP2-SP6) were presented in GOU14. These crucial soft state spectra were all extracted from a single observation, which was 
made on Jan. 6th, 2011. 

Following GOU11 and GOU14, we also corrected for the effective area of the PCA using the Crab Nebula as a calibration 
source (Toor & Seward 1974). We applied several correction factors for RXTE (see Steiner et al. 2010 for specific definitions 
of these correction factors), and the respective values for SP1 and SP2-SP6 are (as described in GOU11 and GOU14): the 
normalization correction factors C's (the ratio of the observed normalization to that of Toor & Seward 1974) of 1.123/1.128; the 


power-law slope correction factors AI pg (the difference between the observed value of T and that of Toor & Seward 1974) of 
0.023/0.022; and the dead time correction factors (the normalization correction factors for detector dead time) of 1.048/1.029. 
For Chandra observation in SP2-SP6, we rebinned the data to achieve a minimum of 200 counts in each new bin. For ASCA 
and RXTE observations, we rebinned the data with at least 25 counts within each bin. 
Finally, we included systematic errors to account for the uncertainties in the instrumental responses: 0.5% for RXTE and and 
1% for ASCA.! 


SPECTRAL ANALYSIS 


We performed our data analysis and spectral fitting using the software package XSPEC version 12.10.17. The fit statistic in use 
is chi-square x?. And the photoelectric absorption cross-sections in use is based on Verner et al. (1996). During the initial fitting 
process, the three key input parameters M, i and D are fixed at their newly-obtained values given in Section 1. 


Our Adopted Model 


Three typical components might be included in the spectrum of a black hole X-ray binary, namely the thermal, power-law and 
reflection components (refer to Figure 2 in GOU11 for a schematic illustration of these different components). 

The model we employed in this paper is the most advanced and physically realistic model used in GOU11 and GOU14. The 
structure of the model, naming all the components that it comprises, with a concise description of each one, are as follows (the 
detailed parameter configuration is the same as GOU11 and GOU14, unless otherwise indicated): 


CRABCOR x CONSTANT « TBABS [SIMPLR ® KERRBB2 + KERRDISK + KERRCONV 8 (IREFLECT & SIMPLC)) 


Thermal component: The thermal radiation is emitted from the accretion disk. KERRBB2 (McClintock et al. 2006) is a fully 
relativistic thin disk model, including the effects of the returning radiation and the limb darkening. This model reads in and uses 
a pair of look-up tables for the spectral hardening factor f, computed using BHSPEC (Davis et al. 2005, Davis & Hubeny 2006) 
and KERRBB (Li et al. 2005), corresponding to two representative values of the viscosity parameter: a = 0.01 and a = 0.1 (here 
we adopt a = 0.1). KERRBB2 only has two principal fit parameters: the dimensionless spin parameter ax, and the mass accretion 
rate M. 

Power-law Component: The thermal seed photons supplied by the accretion disk are inverse-Compton scattered into the power- 
law component by hot electrons in the corona. The convolution model SIMPLR, a modified version of SIMPL (Steiner et al. 20092), 
is an advanced empirical Comptonization model including three main fit parameters: the photon power-law index I’, the scattered 
fraction fsc and the reflection constant parameter x. As we mentioned, we adopt fz. < 25% as our spectrum selection standard. 

Reflection Component: A fraction of the scattered seed photons are emitted downward toward the disk, illuminating the disk 
and generating the reflection component. KERRDISK + KERRCONV® (IREFLECT & SIMPLC)] models the entire reflection com- 
ponent including the Compton hump, the absorption edge and the emission lines. SIMPLC represents the part of the power-law 
component that strikes the disk (Steiner et al. 20092), which is regarded as the input spectrum for the convolution model IRE- 
FLECT (Magdziarz & Zdziarski 1995). Photons emitted near the black hole suffer various relativistic effects, leading to a blurring 
of the observed spectrum. KERRCONV is utilized with the aim of smearing the spectrum (Brenneman & Reynolds 2006). KER- 
RDISK is a model for the iron emission line with the spin allowed to be a free parameter (Brenneman & Reynolds 2006). Here it 
should be noted that the effect on the spin measurements of using different reflection models is almost negligible, as demonstrated 
in section 5.3 of GOU11. 

In addition, we use three multiplicative models: (1) CRABCOR is used for correcting calibration deviations of RXTE. (2) CON- 
STANT is an energy-independent multiplicative factor in order to reconcile calibration discrepancies between different detectors. 
(3) TBABS is an ISM absorption model (Wilms et al. 2000). 


Fit Results 


Our model provides a good fit to all six spectra, with reduced x? values varying from 1.12 to 1.49. The fit parameters are 
shown in Table 2. A representative plot of the unfolded spectrum is given in Figure 1. It is noticed that there is a larger residual 
between 8 keV and 15 keV, which will be discussed in detail later, and should not significantly affect the spin results. 


! To achieve an acceptable fit to the PCA spectra, the standard prescription recommended the addition of a 0.5%-1% systematic uncertainty before fitting the 
spectra, as we did here. However, Garcia et al. (2014) proposed a new empirical calibration model named PCACORR, based on the fitting of hundreds of RXTE 
PCA spectra of the Crab. They found that the quality of the spectra could be significantly improved and required only a 0.196 systematic uncertainty. This was 
tested in García et al. (2014) and García et al. (2016). However, for consistency with the previous measurements in GOUI1 and GOUIA, we did not apply the 
PCACORR tool to the spectra. However, we did one experiment, applying the PCACORR tool to all the PCA spectra and refitting the spectra, and we found the 
results to be almost the same, with negligible changes. 


? XSPEC is available at http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 


As a caveat, in all six cases, the values of a, reached the maximum value of the XSPEC model KERRBB2, 0.9999. This value 
is larger than the physical upper limit of a., which we will discuss in the next section. The preliminary error in a, computed with 
XSPEC and listed in Table 2 is due to counting statistics, while the final adopted error due to the combined uncertainties in D, 
M, and i is acquired via our Monte Carlo analysis. 

Requiring the spectrum to be dominated by accretion disk component for reliable application of the continuum-fitting method, 
we adopt fse < 25% as the selection standard. The values of this parameter in our analysis, 24.7% + 0.5%, 10.6% + 1.6%, 
15.6% + 1.2%, 21.1% + 1.0%, 16.3% + 0.9%, 9.3% + 0.2%, satisfy this standard. 

It is equally essential to select spectra with bolometric Eddington-scaled luminosities | < 0.3 (McClintock et al. 2006, 
Zhu et al. 2012, McClintock et al. 2014), in order to ensure a geometrically thin disk. This parameter is quite stable and low 
in our analysis, ranging from 0.01 to 0.03, which meets this criteria without difficulty. 


Error Analysis 


Previous work on estimating a. had found that the combined observational uncertainties in M, 7 and D dominate the error 
budget (GOU11, McClintock et al. 2014). This part of the error is usually calculated through Monte Carlo simulations (Liu et al. 
2008, Gou et al. 2011, McClintock et al. 2014). In order to do so, for each individual spectrum, the procedure is as follows: 
(1) Assuming that each one of these three parameters is independent and Gaussian distributed, we generate 3000 sets of (M, i, 
D). (2) Then we calculate look-up tables of f for these parameter sets. (3) Finally we fit the spectrum to our adopted model to 
determine a,.. With 3000 simulations run, we obtain the histogram of a... 

After repeating the procedure expressed above for all six spectra, we obtain six histograms of a. The histogram for each 
spectrum is shown in Figure 2, and the summed histogram is plotted in Figure 3. The spin value in SP2 is lowest, therefore we 
report the final result of a, > 0.9985 (3c) based on a conservative estimation from this spectrum. 


DISCUSSION 
Comparison with the Thorne Limit 


In this paper, we have obtained a much more stringent limit on the spin of the black hole in Cygnus X-1 than found from 
previous work: a, > 0.9985 at 3c level of confidence. This very high spin parameter confirms that Cygnus X-1 hosts an extreme 
Kerr black hole. 

However, Thorne (1974) pointed out that the maximum expected spin once the absorption of hot photons from the disk is 
taken into account is 0.998 (hereafter, the Thorne limit), as the absorption will produce an counteracting torque when the spin is 
extremely high. We note that our spin constraint is slightly higher than the Thorne limit. This is mainly because the upper limit 
in our model KERRBB2 is set at 0.9999, and our obtained spin parameter is based purely on the theoretical model. 

It should be noted that although the Thorne limit has generally been adopted by the community as a reasonable guess for an 
astrophysical upper limit, there is still some debate as to the exact value of the upper limit. Sadowski et al. (2011) suggested that 
black holes with super-critical accretion flows could reach an equilibrium spin value up to 0.9994, which exceeds the Thorne 
limit. They showed that the equilibrium spin value depends strongly on the assumed value of the viscosity parameter, and also 
proved that for high accretion rates the impact of captured radiation on spin evolution is negligible. Cygnus X-1 is not believed 
to be in the super-critical regime, and hence these considerations may not apply. However, a number of other factors might affect 
the precise spin measurement, and they are discussed in the following sections. Regardless, the spin of Cygnus X-1 appears to 
be very high, as found previously by both continuum fitting (GOU11, GOU14) and iron line fitting studies (Fabian et al. 2012, 
Tomsick et al. 2014, Parker et al. 2015, Duro et al. 2016, Walton et al. 2016, Tomsick et al. 2018). 


Effect of a Possible Spin-orbit Misalignment 


Figure 4 shows the effect on a, of separately varying the input parameters M, i, and D. a. takes the extreme value of 0.9999 
when all three parameters are fixed on their newly-obtained, best-fitting values. 

In applying the continuum-fitting method, we assume that the binary orbital angular momentum vector is aligned with the 
spin axis of the black hole. Cygnus X-1 is suggested to have formed with little to no natal kick (Mirabel & Rodrigues 2003, 
Rao et al. 2019), so the spin axis of the black hole should still be aligned with the orbital angular momentum vector (as discussed 
by Miller-Jones et al. (2020), who estimated a likely maximum misalignment angle of order 10?). 

However, there is still controversy about this assumption. Some previous work on fitting the reflection component of Cygnus 
X-1 suggested a different inner disk inclination angle compared to the binary orbital inclination, which implies a possible spin- 
orbit misalignment (also suggested in some other black hole binary systems, Fragos et al. 2010, Dong et al. 2020). Tomsick et al. 
(2014) evaluated the lower limit of this misalignment to be 13?. Walton et al. (2016) estimated the magnitude of a disk warp to 
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be about 10° — 15°. However, these fits often require very high iron abundances, and Tomsick et al. (2018) recently found that 
a much lower disk inclination of < 20° could be accommodated if the electron density was free to vary, even suggesting that a 
misalignment might not even be required in this scenario. Thus, while the possibility of a misalignment should be considered, 
this scenario would appear to remain under debate. 

Nonetheless, we investigated this possibility, and as shown in Figure 4(a), the value of a, monotonically decreases with 
increasing i. We find that even given a misalignment angle of 15° (equivalent to i = 42.47°), the spin value is 0.9696 (see Case 1 
in Table 2), which is still very high. 


Effect of the Metallicity 


The spectral hardening factor f in KERRBB2 is pre-calculated with a table model BHSPEC, which is related to the metallicity 
Z. Our fit results are all given for Z = 1 Zo. However, the metallicity of Cygnus X-1 is slightly super-solar (Shimanskii et al. 
2012), as also implied by the fitted super-solar values of Fe abundance in our reflection model IREFLECT. 

Unfortunately, the current model BHSPEC only contains two possible metallicity cases (Z = 1 Zo and Z = 0.1 Zo), and does 
not yet have an option for super-solar abundance. However, if the metallicity increases, the value of the spin increases slightly. 
In addition, substantial changes in metallicity only produce a negligible effect on spin. For example, as stated in Gou et al. 2009, 
the spin of the black hole in LMC X-1 increases from 0.937 + 0.020 to 0.938 + 0.020 if Z increases from 0.1 Zo to 1 Zo. 


Effect of the Viscosity Parameter 


Throughout this work, in order to get a more conservative constraint, we fixed the viscosity parameter a to 0.1, because ax 
increases slightly as a decreases (Gou et al. 2009). However, we reanalyzed our data using a = 0.01, and found that the results 
remain unchanged, as they have already reached their maximum value of 0.9999. 


Effect of a Finite Disk Height 


The model utilized in the continuum-fitting method is the general relativistic Novikov-Thorne model that assumes a razor- 
thin disk. However, a real accretion disk will have finite thickness, which will affect measurements of the black hole spin. 
Kulkarni et al. (2011) evaluated the error introduced by the deviations from the Novikov-Thorne model. They computed disk 
spectra obtained from three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations, and then fitted these 
simulated spectra with KERRBB. By comparing the spin results with the spin used in the GRMHD simulations, they concluded 
that the model errors become obvious only at low spins and high inclinations. In the case of ? = 30° and a, = 0.98 (close to the 
values estimated for Cygnus X-1), the spin calculated by fitting the simulated spectra was 0.985 + 0.001 and the value of A 
ax is 0.005. Moreover, they might have overestimated the error due to the deviation from the Novikov-Thorne model because 
their simulation results are for a thicker disk with a bolometric Eddington-scaled luminosity / = 0.4—0.7, while the bolometric 
Eddington-scaled luminosity of Cygnus X-1 in this paper is only 0.02 (as stated in Section 6 of GOU11). 


The Residual Structure between 8—15 keV 


We notice that there is a larger residual between 8 keV and 15 keV, which is not fitted with the current model. However, it 
should be emphasized that the spin parameter is determined from the disk thermal component (which dominates the spectra, 
peaks around 1.5 keV, and is concentrated below 6 keV; see Figure 3 of GOUI1 for reference), and thus the residual between 8 
and 15 keV should not significantly affect the spin results. 

Previous work has proposed many theories to explain similar residual structures. Fabian et al. (2020) found an 7-9 keV 
excess emission component during the soft state in spectra of MAXI J1820+070. They explained this additional component by 
including another blackbody whose temperature is slightly higher than that of the disk itself. Given that the higher temperature, 
the additional emission is linked to the innermost region of the disk and with the start of the plunging region. In this region, the 
density is still high enough to heat the gas, possibly emitting blackbody radiation. We utilized an addition model BBODY to 
describe this component, we find that the fit is poor (> 2), and the spin could not be determined certainly. 

Another possibility is that the observed residual structure is an additional Comptonization component. In our fits, we only 
utilized Comptonization models which have a pure non-thermal power-law electron distribution. However, some previous work 
of the fast time variability seen in the high/soft state prefer a stratified Comptonization spectrum (Grinberg et al. 2014). It 
is likely that there exist hybrid thermal/non-thermal electron distributions, resulting in an additional low temperature thermal 
comptonization component (?, Zdziarski et al. 2002, Kawano et al. 2017). This component then produce a higher energy power- 
law tail, which would also lead us to overestimate the black hole spin. For instance, as mentioned in Kawano et al. (2017), the 


spin of the black hole in Cygnus X-1 is 0.95 + 0.01 for a simple non-thermal Comptonization; and decreases to 0.80 ure 
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for hybrid thermal/non-thermal Comptonization. A similar residual structure between 5—10 keV was found in the spectra of 
the black hole candidate MAXI J1828-249 in a possible intermediate state (Oda et al. 2019). They modeled this excess with a 
Comptonization of disk photons by thermal electrons with a relatively low temperature and inferred it could be generated in a 
region where the disk inner edge intrudes into the hot inner flow, or a Comptonizing region above the disk surface. We used an 
addition model COMPTT to describe the Comptonization of soft photons, and reflected both two Comptonization components. 
The model is reconstructed as follows: 


CRABCOR « CONSTANT « TBABS [KERRBB2 + SIMPL & COMPTT 
+KERRDISK + KERRCONV @ (IREFLECT & SIMPL & COMPTT)] 


The fit parameters are listed in Table 3. For SP1, the fit is better than former model with reduced x? reducing from 1.49 to 
0.94. The spin parameter a, decreases from 0.9999 to 0.9947 + 0.0014. The photon power-law index decreases from 2.39 + 0.01 
to 1.99 + 0.01. The scattered fraction fsc is larger ranging from 0.247 + 0.005 to 0.399 + 0.007. For SP2-SP6, the values of ax 
are unaffected, which are still at 0.9999. However, the errors quoted here (estimated in XSPEC) are obviously larger. Besides, 
for these five spectra, this model can't constrain the parameters of COMPTT well probably because the component is weak. The 
scattered fraction fsc are somehow increasing from ~ 0.15 to ~ 0.40. For SP5, the fit is slightly worse with reduced x? increasing 
from 1.25 to 1.27. While for SP2, SP3, SP4 and SP6, the fit is better with reduced x? ranging from 1.12, 1.15, 1.32, 1.24 to 0.93, 
1.02, 0.94, 1.12, respectively. We also try the same setup of Kawano et al. (2017), fixing To, kT and 7 at 0.44 keV, 3.4 keV and 
1.66, respectively. Details are shown in Case 1 of Table 3. The spin parameter a. is still pegged at 0.9999. 

Whether using a hybrid thermal/non-thermal Comptonization model or a single non-thermal Comptonization model, we all get 
an extremely high spin. So we can conclude that including COMPTT is not having the same effect as in Kawano et al. (2017). 
However, it is interesting to note that the reduced x? decreases significantly (e.g., for SP1 it goes to 0.94 from 1.49), which means 
that adding this Comptonization model is having a large impact on the high energy Compton component rather than on the disk 
which is the foundation of our spin measurements. 


Comparison with the Reflection Measurements 


Another widely-used technique of measuring the spin is the iron-line profile fitting method (Fabianetal. 1989, 
Reynolds & Nowak 2003). As mentioned above, some portion of the power-law component is reflected by the disk, generat- 
ing the reflection spectrum whose main prominent features include the fluorescent emission lines (in particularly, the iron Ko 
line between 6.4 - 6.97 keV, depending on ionization state) and the Compton hump at 20 - 30 keV (Ross et al. 1999). The in- 
trinsically narrow iron Ko line is then broadened and smeared by the combined effects of Doppler shift, special relativity, and 
general relativity (Fabian et al. 2000), and by fitting the iron-line profile, one can deduce the spin parameter. 

The spin results from the continuum-fitting method both in the previous work and in this work are consistent with the spins 
derived from the reflection method (Duro et al. 2011, Fabian et al. 2012, Tomsick et al. 2014, Parker et al. 2015, Duro et al. 2016, 
Walton et al. 2016, Tomsick et al. 2018), although Miller et al. (2002) found a near-zero spin of Cygnus X-1, which is mainly 
because the lack of adequate energy coverage (< 10 keV) made it difficult to subtract the underlying continuum properly. 

It should be noted that the highly-sensitive NuSTAR, launched in 2012, provides unprecedentedly precise measurements of 
the iron line region and reflection component. In the soft state, Tomsick et al. (2014) reported a rapid spin of a. > 0.83 and 
Walton et al. (2016) estimated a higher spin of 0.93 < a, < 0.96. Parker et al. (2015) analyzed observations taken in the hard 
state and found a high spin of a, > 0.97. In the intermediate state, Tomsick et al. (2018) restricted the spin to be a, > 0.93 for 
the high-density reflection model. 

In addition, Duro et al. (2016) analyzed a spectra of four simultaneous hard intermediate state observations, and found a spin 
of ~ 0.9. They also suggested that the different inclination angles inferred in the different states (i.e., the reported inclination 
angle of ~ 30° in the hard intermediate state, as compared to the higher inclinations found from the reflection measurements in 
the soft state) could imply a warp in the inner region of the disk. However, as discussed by Tomsick et al. (2018), these results 
can be somewhat degenerate with the iron abundance and with the disk electron density. 


Spin Comparison with Other Systems 


We found that the spin of Cygnus X-1 is constrained to be a. > 0.9985, which is the highest among all the stellar-mass black 
hole systems with a measured spin parameter. Although it may be less extreme due to some systematics, this value would be 
still very high. If we look through the spin parameters measured for the stellar-mass black hole systems using the continuum- 
fitting method, it shows that almost all the black holes in high-mass X-ray binaries have a relatively high spin (Cygnus X-1: 
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ax > 0.9985, this work; LMC X-1: a, = 0.92+9:05, Gou et al. 2009; M33 X-7: a, = 0.84 + 0.05, Liu et al. 2010; IC10 X-1: 
ax Z 0.8, Steiner et al. 2016). However, the low mass X-ray binaries have a wide range of spins between zero and unity, which 
may imply different evolutionary paths for these systems (Bambi 2019, Reynolds 2019). 

It is noted that, since the first detection of GW150914 in 2015, LIGO and Virgo have made 3 runs, and also made 50 detections 
of gravitational wave events with majority of stellar mass binary black hole systems. The results from the first 2 runs have been 
catalogued in Abbott et al. (2019), which shows the effective spin Xe *of the binary system before merger is quite small with the 
distribution’s mean u ~ 0, and with a narrow distribution, which is confirmed by Miller et al. (2020). Miller et al. (2020) further 
explored what these ensemble properties imply about the spins of individual binary black hole mergers, re-analyzing existing 
gravitational-wave events with a population-informed prior on their effective spin, and found, under this analysis, the binary 
black hole GW170729, which previously excluded ye = 0, is now consistent with zero effective spin at ~ 10% credibility. 

As we know, the companion star of Cygnus X-1 is around 40 Mo, and based on our understanding to the stellar evolution, 
it may eventually produce a second black hole of the system, forming a binary black hole system in the future. In addition, the 
current analysis to the system shows that the black hole spin-orbit misalignment is small, within 10 degree (Miller-Jones et al. 
2020). If the black hole of Cygnus X-1 and the companion star were formed from the same cloud initially, then the companion 
star is usually aligned with the orbital plane. In this case, we should expect the effective spin of the system is quite large, unless 
there is a large asymmetric kick during the black hole formation of the companion star, and it will be quite different from the 
spin distribution obtained from the GW results. This may imply that Cygnus X-1 system is likely to follow a different formation 
channel (primordial one) compared to the GW black hole system, although it is not clear whether the GW results preferentially 
support a dynamical or field binary origin (or even other possibilities) for stellar mass binary black holes detected by GWs. 
Furthermore, Miller-Jones et al. (2020) shows that, given the current orbital separation, we do not expect Cygnus X-1 or systems 
with similar properties to yield binary black hole mergers in the age of the Universe. So we may haven't detect any GW signal 
from primordial binary black hole systems yet. 

As to the origin of the high spin in these high-mass systems, it is usually concluded that the spin is natal (Podsiadlowski et al. 
2003, Axelsson et al. 2011, Qin et al. 2019), at least for those high-mass X-ray binaries, mainly because in the whole lifetime 
of a stellar-mass black hole, it is centainly impossible to spin up from zero to a highly-spinning black hole via standard (as 
opposed to hypercritical) accretion (Wong et al. 2012). One possible evolutionary channel to explain the high spin is that the 
BH progenitor undergoes mass transfers while still in the main sequence (see Case-A MT in Qin et al. 2019). Hence the core 
of the BH progenitor contains enough angular momentum to form a highly rotating black hole (many observational features 
predicted by this path, such as unusual surface abundances of the secondary, have been found in Cygnus X-1, as mentioned in 
Miller-Jones et al. 2020). 


CONCLUSION 


We re-analyze six archival spectra of the black hole X-ray binary Cygnus X-1, which were originally presented in GOU11 and 
GOU14, to constrain the black hole spin. We adopt up-to-date values of three key dynamical parameters M, i, and D. These data 
rigorously satisfy the selection criterion fse < 25%. The key model we utilized is a fully relativistic thin disk model KERRBB2. 
Monte Carlo simulations are performed in order to estimate the error in a, due to the combined uncertainties of M, i, and D. 
Ultimately, we arrive at our final result, a, > 0.9985 at the 3c level of confidence, which confirms the extreme spin of the black 
hole in Cygnus X-1. 

It should be noted that the spin value we presented above is the theoretical one purely based on the model itself, and the hard 
limit is set at 0.9999. In practice, the maximum spin value could be slightly lower than 0.9999. Thorne (1974) estimated the 
evolution of a black hole, and pointed out that the maximum spin of a black hole can only reach 0.998 because the radiation from 
the accretion disk will produce a counteracting torque on the black hole once the black hole reaches the maximum spin. Given 
the Thorne limit, we conclude that the black hole in Cygnus X-1 has an extreme spin in any event. 
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"Xe = T Un ami) 0 
angular momentum. 


, where m; and d; indicate masses and dimensionless spins of each component, Lx is the unit vector along the direction of the orbital 
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Number 


Energy Band (keV) 


ObsID 


Mission 


Table 1. Log of the Observations? 


Detector 


UT Texp (S) —I(Crb) ¢ 
SP1 10408000 & P10412 ASCA & RXTE GIS & PCA 0.7-8.0 & 2.5-45.0 1996-05-30 06:43:16 & 07:51:29 2547 & 2240 0.80 — 0.74 
SP2 12472 & P96378 Chandra & RXTE | HETG(CC) & PCA 0.8-8.0 & 2.9-50 2011-01-06 14:06:40-14:35:44. — 455 & 1744 0.32 0.32 
SP3 12472 & P96378 Chandra & RXTE | HETG(CC) & PCA 0.8-8.0 & 2.9-50 2011-01-06 15:44:16-16:09:52 398 & 1536 — 0.61 — 0.33 
SP4 12472 & P96378 Chandra & RXTE | HETG(CC) & PCA 0.8-8.0 & 2.9-50 2011-01-06 17:15:28-17:43:44  455& 1744 0.57. 0.35 
SP5 12472 & P96378 Chandra & RXTE | HETG(CC) & PCA 0.8-8.0 & 2.9-50 2011-01-06 18:19:44-19:17:52 455 & 1744 038 — 0.36 
SP6 12472 & P96378 Chandra & RXTE | HETG(CC) & PCA 0.8-8.0 & 2.9-50 2011-01-06 19:53:36-20:50:08 455 & 1744 038 0.37 


NOTE—“ The log of the observations analyzed in this work (adapted from observing logs in GOU11, GOU14). In columns 2-9, we show 
the following information, respectively: observation identifier; names of the observational missions; names of the detectors requested for the 
exposure; energy band used for the spectral fit; the start and end times of the observations; the effective exposure times; the intensity of the 
source; The orbital phase of the binary system, which is defined (at the midpoint of the observation) relative to the time of superior conjunction 
of the O-star (black hole beyond star), occurring on heliocentric Julian Day 2441163.529 (Orosz et al. 2011). 


Table 2. Best-fit Results for Our Adopted Model" 


Spectral Number 


Number Model Parameter? Case 1 
SPI SP2 SP3 SP4 SP5 SP6 
1 TBABS Nu 0.74 + 0.01 0.82 + 0.02 0.83 + 0.03 0.80 + 0.02 0.79 + 0.02 0.81 +0.01 0.72 +0.01 
2 SIMPLR r 2.39 + 0.01 2.38 + 0.04 2.46 + 0.03 2.69 + 0.02 2.62 + 0.03 2.50 + 0.01 2.45 +0.01 
3 SIMPLR p 0.247 + 0.005 0.106 + 0.016 0.156 + 0.012 0.211 + 0.010 0.163 + 0.009 0.093 + 0.002 0.092 + 0.002 
4 KERRBB2 as 0.9999 0-0009 0.9999 0-0099 0.9999: 0:0999 0.99997 9:6090 0.9999:0-0000 0.9999 00015 0.9696 + 0.0034 
5 KERRBB2 M 0.17 + 0.02 0.20 + 0.03 0.19 + 0.04 0.19 + 0.03 0.18 + 0.03 0.18 + 0.01 0.23 + 0.01 
6 KERRDISK E, 6.57 + 0.06 6.74 + 0.08 6.55 + 0.07 6.40 + 0.13 6.55 + 0.11 6.70 + 0.08 6.46 + 0.06 
j KERRDISK q 2.89 + 0.08 3.00 + 0.09 2.73 + 0.11 1.64 + 0.40 2.66 + 0.17 2.95 + 0.06 345 + 0.10 
8 KERRDISK NL 0.009 + 0.001 0.024 + 0.003 0.029 + 0.003 0.012 + 0.002 0.011 + 0.002 0.014 + 0.001 0.026 + 0.001 
9 KERRDISK EW 0.11 0.36 0.29 0.14 021 0.29 0.61 
10 IREFLECT Xre 6.3 + 0.3 6.01.1 3.9 +0.5 3.5 + 0.3 5.9 + 0.6 4.1 +0.6 3.5 + 0.6 
11 IREFLECT € 62 3-10 335 + 186 137 45 2147 57418 381 +51 412 + 64 
12 CONSTANT Nasca 0.955 + 0.003 - - - - - - 
13 CONSTANT Nec = 0.688 + 0.009 0.575 + 0.007 0.704 + 0.007 0.713 + 0.005 0.682 + 0.006 0.690 + 0.011 
14 CONSTANT Nusa < 0.863 + 0.011 0.889 + 0.012 0.887 + 0.009 0.923 + 0.007 0.889 + 0.008 0.903 + 0.014 
15 Reduced x? (x?/degree) 1.49 (678.09/454) 1.12 (705.15/627) 1.15 (656.82/572) 1.32 (823.61/624) 1.25 (929.81/744) 1.24 (1241.40/998) — 1.09 (1087.79/998) 
16 f 1.6 16 1.6 16 1.6 16 15 
17 l 0.02 0.03 0.03 0.02 0.02 0.02 0.02 


NOTE—" Names of the different model components are listed in column 2. In column 3, the parameters from top to bottom are: (1) equiv- 
alent hydrogen column density Ny in units of 107? atoms cm ?; (2) the photon power-law index T; (3) the scattered fraction fsc; (4) the 
dimensionless spin parameter a, (the upper limit is set to 0.9999); (5) the mass accretion rate M in units of 1018 g s7}; (6) the rest frame 
central line energy £r, in units of keV; (7) the emissivity index q; (8) the line flux Nz, in units of photons cm? s7}; (9) the equivalent width 
EW of the line in units of keV; (10) the iron abundance Xre, relative to the solar abundance; (11) the disk ionization parameter £ in units of 
ergcm s! (12) the detector normalization constant Nasca of ASCA relative to RXTE PCU2; (13) the detector normalization constant NMEG 
of Chandra MEG relative to RXTE PCU2; (14) the detector normalization constant Nueg of Chandra HEG relative to RXTE PCU2; (15) the 
reduced chi-squared value, total chi-squared value and number degrees of freedom; (16) the spectral hardening factor f; (17) the bolometric, 
Eddington-scaled X-ray luminosity of the thermal disk component, | = L/Lzaa, where L = L(ax, M ) and Lgaa = 2.8 x 108° erg s+ for M = 
21.2 MO. The uncertainties are quoted at the 90% confidence level. 

"Case 1: fit to SP6 with the inclination angle i fixed at 42.47°. The parameter configuration is exactly the same as that in column 5-9. 


Number Model 


1 TBABS 
2 KERRBB2 
3 KERRBB2 
4 SIMPLR 
5 SIMPLR 
6 COMPTT 
T COMPTT 
8 COMPTT 
9 COMPTT 
10 KERRDISK 
11 KERRDISK 
12 KERRDISK 
13 KERRDISK 
14 IREFLECT 
15 IREFLECT 
16 CONSTANT 


17 CONSTANT 
18 CONSTANT 


Table 3. Best-fit Results for a Hybrid Thermal/Non-thermal Comptonization Model 


Parameter 


Xre 
é 
Nasca 
Nec 


Nuec 


Reduced x? (x?/degree) 


f 
L 


Spectral Number 


SP1 


0.730.01 


SP2 


0.73 + 0.05 


0.9947 + 0.0014 0.9999 +0-0000 


0.18 + 0.01 
1.99 + 0.01 
0.399 + 0.007 
0.453-0.06 
5.7742.51 
2.08+0.74 
1.0840.48 
6.59+0.04 
2.74+0.03 
0.019+0.001 
0.21 
2.2+0.1 
19811 
1.012+0.003 


0.14 + 0.06 
2.40 + 0.04 
0.447 + 0.440 


0.45 + 0.44 
4 


4 


+37.60 
0.81*37.6 


2.50* 5509 
6.56 + 0.07 
2.51 £0.19 
0.020 + 0.004 
0.29 
5312 
107 + 48 
0.824 + 0.010 
1.037 + 0.012 


0.94 (430.62/451) 0.93 (581.00/624) 


1.6 
0.02 


1.6 
0.02 


SP3 


0.71 + 0.06 
0.9999 02105 
0.22 + 0.02 
2.46 + 0.02 
0.431 + 0.359 
0.40 + 0.23 
4.00* 115,7? 
0.63 5 
3.64* 15,72 
6.51 + 0.08 
2.50 + 0.21 
0.023 + 0.004 
0.25 
3.8 + 0.6 
80 +31 
0.664 + 0.008 
1.029 + 0.012 
1.02 (580.58/572) 
L5 
0.01 


SP4 


0.65 + 0.06 
0.999970 1098 
0.25 + 0.03 
2.53 + 0.02 
0.393 + 0.198 
0.41 + 0.36 


6.49 + 0.08 
2.52 + 0.22 


0.020 + 0.004 
0.23 
4.5 x: 0.8 
92 +37 
0.813 + 0.009 
1.026 + 0.012 
0.94 (582.79/621) 
1.6 
0.02 


SP5 


0.71 + 0.03 
0.9999 0.0009 
0.13 £0.10 
243 + 0.05 
0.361 + 0.234 


0.39 + 0.16 
4.067388 


2.158537 
6.57 + 0.07 
2.53 + 0.19 
0.014 + 0.003 
0.29 
414 10 
111 +48 
0.803 + 0.007 
1.039 + 0.009 
1.27 (941.07/741) 
1.6 
0.02 


SP6 


0.75 + 0.04 
0.9999 0:0099 
0.29 + 0.05 
2.52 + 0.02 
0.361 + 0.311 


0.42 +0.26 
395% sas. 
0.43 3 
2417522 
6.55 + 0.08 
2.57 £0.19 
0.013 + 0.002 
0.26 
3.9 +08 
128 +53 
0.800 + 0.008 
1.045 + 0.010 
1.12 (1114.26/994) 
16 
0.02 


Case 1 


0.76 + 0.02 
0.99995 00009 
0.16 + 0.03 
2.39 + 0.02 
0.492 + 0.176 
0.44 (fixed) 
3.4(fixed) 
1.66(fixed) 
2.12+ 0.07 
6.54 + 0.07 
2.38 + 0.16 
0.019 + 0.002 
0.26 


0.820 + 0.009 
1.031 + 0.011 
0.93 (583.43/627) 
1.6 
0.02 
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NOTE—^ Names of different parameters in model COMPTT are (1) the input soft photon Wien temperature To in units of keV; (2) the plasma 
temperature kT in units of keV (the lower limit is set to 0.1 keV); (3) the plasma optical depth 7 (the lower limit is set to 0.01); (4) the 
normalization norm (the lower limit is set to 0). 


"Case 1: fit to SP2 with To fixed at 0.44 keV, kT fixed at 3.4 keV and 7 fixed at 1.66. The parameter configuration is exactly the same as that 


in column 5-9. 


Photons cm? s- keV-! 


(data-modely'error 


Unfolded Spectrum 
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"^g 


Energy (keV) 


20 50 


Figure 1. A representative best-fit X-ray spectrum (SP3). Left: A plot of the unfolded spectrum and the model. Right: A plot of the spectrum 
with the folded model. Detailed fit parameters are listed in Table 2. The data obtained by RXTE is plotted in red. For the Chandra observation we 
used, the data collected by the High Energy Grating (HEG) and Medium Energy Grating (MEG) components of the High-Energy Transmission 
Grating (HETG) are plotted as green and black, respectively. The dotted lines represent the models we adopted. 
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a->0.9995 


(SP6) 
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N 
] 


id 
[n] 
1 


Probability 


0.34 


I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
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0.2 4 | 
I 
I 
I 
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0.9988 1.0000 


0 T 
0.9977 0.9988 


Figure 2. Histograms of a. calculated via the Monte Carlo analysis for 3000 sets of parameters per spectrum. The three dashed lines from left 
to right imply the 99.7% (3c, red), 95.4% (2c, yellow), and 68.396 (1e, green) error, respectively. The respective 99.7% confidence lower limit 
on a. is indicated in each panel. 


ax>0.9987 


0.9984 0.9986 0.9988 0.9990 0.9992 0.9994 0.9996 0.9998 1.0000 
ax 


Figure 3. Summed histogram of a. for our six spectra, including 18000 data points. The dashed lines and listed 99.7% confidence lower limit 
in this figure are as those shown in Figure 2. 
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Figure 4. Effects on a. of varying one of the three input parameters while fixing the other two at their updated values. (a) Spin vs. inclination 
angle i; (b) spin vs. black hole mass M; (c) spin vs. source distance D. The red dashed lines represent the baseline values i = 27.51 °, M = 
21.2 and D z 2.22 kpc. These plots were made using the fit results of SP6. 


